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Abstract 

Communication between cells is a ubiquitous feature of cell populations and is frequently realized by secretion and 
detection of signaling molecules. Direct visualization of the resulting complex gradients between secreting and receiving 
cells is often impossible due to the small size of diffusing molecules and because such visualization requires experimental 
perturbations such as attachment of fluorescent markers, which can change diffusion properties. We designed a method to 
estimate such extracellular concentration profiles in vivo by using spatiotemporal mathematical models derived from 
microscopic analysis. This method is applied to populations of thousands of haploid yeast cells during mating in order to 
quantify the extracellular distributions of the pheromone a-factor and the activity of the aspartyl protease Barl. We 
demonstrate that Barl limits the range of the extracellular pheromone signal and is critical in establishing a-factor 
concentration gradients, which is crucial for effective mating. Moreover, haploid populations of wild type yeast cells, but not 
BAR! deletion strains, create a pheromone pattern in which cells differentially grow and mate, with low pheromone regions 
where cells continue to bud and regions with higher pheromone levels and gradients where cells conjugate to form 
diploids. However, this effect seems to be exclusive to high-density cultures. Our results show a new role of Barl protease 
regulating the pheromone distribution within larger populations and not only locally inside an ascus or among few cells. As 
a consequence, wild type populations have not only higher mating efficiency, but also higher growth rates than mixed 
MATa barlA/MATa cultures. We provide an explanation of how a rapidly diffusing molecule can be exploited by cells to 
provide spatial information that divides the population into different transcriptional programs and phenotypes. 
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Introduction 

Cells communicate with each other by detecting and responding 
to external cues and stimuli. In cellular systems one can find 
several examples where cells coordinate growth in specific regions 
by sensing and responding to small differences in the concentra- 
tions of substances that provide positional information [1,2]. Such 
signaling among cells is a common feature of cell populations as 
well as multicellular organisms, where cells often operate close to 
the physical limit of gradient or concentration detection [3-5] . 

The pheromone response of budding yeast (Saccharomyces 
cerevisiae) is an example for a eukaryotic cell communication 
system. Yeast cells occur either in the haploid forms MATa and 
MATa, or as MATa./a diploid. Haploid and diploid cells are both 
able to replicate vegetatively. Mating of two haploid cells with 
opposite mating types yields diploid cells, while haploid cells are 
formed through spore formation in meiosis [6,7]. 

Mating is initiated by the secretion of mating type-specific 
pheromones, called a-factor and a-factor, which are sensed by 



haploid cells of the opposite mating type and trigger the mating 
response [8,9] . During the mating response, yeast cells arrest their 
cell cycle in Gl phase and elongate in the direction of the 
pheromone signal by forming directed mating projections called 
"shmoos" [6,10]. Yeast can sense pheromone gradients as well as 
absolute concentration levels. Nevertheless, yeast cells are not 
capable of chemotaxis and, thus, mating requires the haploid cell 
types to signal their location particularly to nearby potential 
mating partners. One way for MATa cells to regulate the 
extracellular a-factor is the secretion of Barl, an aspartyl protease, 
which degrades a-factor [11,12]. This leads to the paradox that 
MATa. cells degrade the signal they need to receive. Theoretical 
investigations hypothesized that the major role of Barl is to 
sharpen the pheromone gradient [13,14]. 

However, direct visualization of the resulting spatial a-factor 
concentration profiles between secreting and receiving cells is 
impossible due to the small size of diffusing molecules and because 
such visualization requires experimental perturbations such as 
attachment of fluorescent markers, which change diffusion 
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Author Summary 

Haploid budding yeast cells cannot actively move to find a 
mating partner, like some flagellated bacteria do. Instead 
they must grow a so-called shmoo - a mating projection - 
precisely into the direction of a potential partner. They 
communicate with each other by releasing pheromones 
into their environment, which are sensed by cells of the 
opposite mating type. This serves the localization of 
nearby cells and initiates growth arrest and mating. 
Paradoxically, yeast cells also secrete the protease Barl 
that destroys pheromones. To visualize the resulting 
pheromone distribution and understand the effect on 
mating efficiency, we combined fluorescence imaging and 
mathematical modeling. We observed that the controlled 
destruction of pheromones by the yeast cells is beneficial 
to communication since it causes relatively higher pher- 
omone concentrations in areas where cells are dense and 
vanishing pheromone concentrations elsewhere. This 
allows the population to maintain two different cellular 
behaviors at the same time, i.e. mating and continued 
growth, a behavior which disappears when we genetically 
delete the gene for the pheromone-destroying protein. 



properties and activities. Therefore, reaction-diffusion (RD) 
models have been used to simulate the pheromone distribution 
on the basis of physical properties of the molecules [13-16], but 
neither the validity of the model predictions nor the effect of the 
pheromone distribution have been tested experimentally. Except 
of the publication of Jin etal. [15], where it has been suggested that 
Barl promotes avoidance of the same mating type and accurate 
gradient detection. Using a microfluidic device it had been shown 
how MATa cells avoid each other when exposed to an artificial 
unidirectional gradient, which was reproduced quite vividly by 
simulations of an RD model. However, the assumptions and 
choice of parameters were in contrast to other works [14,16] 
(compare Text SI). In general, experimental validation of RD 
models is complicated, especially on the molecular level [17]. 

Moreover, recent theoretical findings support the theory, that 
secretion of Barl in the extracellular medium does not help to 
align gradients in the direction of the opposing mating type [16]. 
In summary, the role of Barl is still controversially discussed. Also, 
none of the models proposed so far has investigated interactions of 
more than a few cells or the four haploid spores in an ascus, even 
though mating occurs not only inside the ascus, but also in a cell 
population which was shown by new findings where a remarkably 
high outcrossing rate from asci was reported [18]. This indicates 
that mating yeast cells interact with quite a number of potential 
mating partners in a natural environment. Furthermore, a recent 
study has shown the potential of simple secrete and sense motifs to 
exhibit surprising effects on the population level [19]. 

Therefore, we designed a method to identify the most likely a- 
factor distribution within mixed haploid yeast populations of 
thousands of cells directly from confocal microscopic images with 
fluorescently tagged marker proteins. Here, an RD model is used 
to simulate interactions of a few hundred cells at the same time. In 
MATa. cells, the protein Fusl, which is strongly expressed upon 
pheromone stimulation, is tagged with GFP to record pheromone 
pathway activation and serves as a proxy for Barl induction. 
Therefore, the experimentally observed pheromone activation 
level of each MATa cell is integrated into the model and compared 
to the experiments. MATa cells are modified to express mCherry 
from the TDH3 promoter to indicate their mating type and 
location. In a unique way we coupled physical RD models with 



experimental imaging in order to quantify the spatial distribution 
of extracellular a-factor. The use of simple marker constructs, that 
altered neither a-factor nor Barl, served for minimal interference 
with the biological system. We used this approach to direcdy 
estimate the influence of Bar 1 on the distribution of a-factor in a 
mixed yeast population and suggest a novel function of Barl to 
enable the coordination of mating and growth in a yeast 
population in vivo. 

Results 

Extracellular a-factor concentrations were quantified 
along with Barl activity directly from in vivo confocal 
images 

We combined image analysis with spatiotemporal mathematical 
modeling to determine spatial concentration distributions of Bar 1 
and of a-factor. Figure 1 introduces the concept of the approach: 

(i) Take images, detect cell location and mating type, and 
quantify pheromone stimulation of MATa. cells with a 
fluorescent marker. 

(ii) Use a mathematical model based on real cell location and 
activation to calculate the distribution of Bar 1 and a-factor 
in the extracellular space. 

(iii) Predict effects on population growth and mating efficiency, 
confirmed by further experiments. 

In order to obtain in vivo conditions during microscopy, which 
were suited for the described methodology, confocal microscopic 
images were taken from synchronized haploid cells or from equally 
mixed haploid MATa,/ MATa cell populations. The cell culture 
samples were spun down with low g-force on glass bottom dishes in 
order to have immobile cells without using concanavalin A coated 
dishes (which we found to alter the population response). This 
protocol essentially yielded sedimented cells in the same way as 
they would be present in any laboratory or naturally occurring 
non-agitated medium, with the difference that sedimentation here 
was achieved under controlled conditions. 

Location, shape, and mating type of the cells were extracted 
from out-of-focus images in the brightfield and mCherry channel 
(Figure SI in Text SI). The cells' spatial arrangement on the 
images was transferred into locally refined triangular meshes for 
the model (Figure S2 in Text SI) and used to calculate the 
extracellular spatial distributions of Bar 1 and a-factor. 

In mathematical terms the problem (ii) is described as a pure 
extracellular reaction-diffusion process for a-factor and Barl with 
distinct boundary conditions (Figure S3 in Text SI). We 
formulated an RD model by the following set of partial differential 
equations: 

Botix t) 

— =D x Acc(5t,t)—B0c,t)-a:(5c,t), 
at 

8B ^=D B AB m . 

This system covers the extracellular dynamics of a-factor 
concentration, cc(x,t), and the activity of Barl protease, B(x,t), 
over time t and position in space, x. The equations quantify two 
types of processes: (1) diffusion of both a and B, where D a and Db 
are diffusion constants and A denotes the Laplacian, and (2) 
degradation of oc by B. Boundary conditions at the cell surfaces 
define secretion of a-factor by MATa cells and induced secretion 
of Barl by MATa cells. 
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Figure 1 . The combined experimental and computational methodology used to derive the extracellular distribution of a-factor. (A) 

Cells were synchronized by elutriation, fixed by forced sedimentation, and microscopically analyzed. (B) Fus1-GFP was measured in MATa wild type 
and bar1Ace\\s after stimulation with various levels of a-factor (displayed here: a-factor concentration of 500 nM) and used for calibration curves in 
(D) and (E). (C) Overlays of the bright field, mCherry and GFP channels allowed recording the positions of MATa (in red) and MATa cells for further 
mathematical analysis. (D) Calibration of the response of barl/t cells to given amounts of a-factor allowed to estimate, how much a-factor would 
reach the cell without degradation by Bar1. (E) Calibration of the response of wild type MATa cells enabled calculation of Bar1 secretion in response 
to a-factor. (F) The information obtained in panels (B) to (E) was used to calculate the unknown distributions of Bar1 and a-factor in the space 
between the cells (a-factor distribution is displayed). 
doi:10.1371/journal.pcbi.1003690.g001 



The diffusive flux of a-factor on the cell surface and exterior 
system boundaries is given by: 

{-J^ ec , on MATa cells 
0 , on MATa cells 
0 , otherwise. 

Here, J?" c is a constant specifying the secretion of a-factor. The 
vector ft points towards the cell interior and, therefore, the 
diffusion flux D c ,ya(x,t)-n at the boundary takes a negative value 
for the secretion of molecules and a positive sign for the 
absorption of molecules. The induction of Bar 1 -activity at each 
MATa. cell is calculated from the average a-factor concentration 
at the surface of this cell, which is represented by a Hill-curve 
(compare Figure 1 E, for details see Figures S4 and S5 in Text 
SI): 



Jg,0 -*»+*! , mXJ ~J )] r- 

\Ii{a{x,t-T)} B +EC H ) 



The expression Ij[tx(x,t — t)] is the average concentration of a- 
factor at the z'-th MATa cell at time t — T, which promotes the 
secretion of Barl. We assume a zero Barl flux on the MATa 
cells and the system exterior: 

D B VB(x,t)n = 0. 

The remaining exercise was to identify the parameter values of 
the RD model. Diffusion constants for the proteins were directly 
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calculated from protein properties (size and density). Thus, only 
two parameters had to be identified by parameter estimation: 
the activity of Barl and the secretion rate of a-factor, as 
described below. 

To determine the activity of Barl, we used a three-step 
procedure. First, we performed an initial calibration to quantify 
the a-factor concentration perceived by individual MATa cells. 
We used MATa. cells carrying a BAR1 deletion (barl A) and the 
pheromone response marker Fus 1 -GFP and synchronized them in 
Gl phase, where they are responsive to pheromone. Their 
response to varying concentrations of a-factor (artificially added, 
in the absence of MATa cells) was quantified as Fusl-GFP 
fluorescence intensity on microscopic images. Since barl A strains 
do not secrete the protease, the local a-factor concentration was 
equal to the applied concentration. Fluorescence intensity of Fusl- 
GFP in correlation with a-factor concentration was recorded as a 
calibration curve (see Figure S4 in Text SI). The calibration curve 
was then applied to mixed haploid cultures to determine the 
perceived a-factor concentration for each MATa. cell on the image 
(see Figure 1). In mathematical terms, the calibration yielded the 
boundary value of a-factor concentration at the surface of a MATa 
cell and a functional relation between this a-factor concentration 
and the induced Barl expression. 

Second, the steady state activity of Barl was quantified by 
stimulating wild type MATa cells carrying the Fusl-GFP marker in 
G 1 phase with given concentrations of a-factor. Since there were 
no MATa cells, the a-factor secretion rate was equal to zero at this 
stage. For the model, these data were used in a mathematical 
optimization to quantify the activity of Barl for wild type cells. In 
order to verify that Fusl can in fact be used as a proxy for Barl, 
we used a strain expressing qVenus fluorophore under the control 
of the Barl promoter (Barlpr-qVenus, [20]) while maintaining 
wild type Barl activity (see Figure S5 in Text SI and Materials 
and Methods). The induction occurred with the same kinetics and 
Hill-coefficients as the Fusl induction, but with a slight delay of 
the Barl expression verifying our used induction kinetics for Barl. 
The Barlpr-qVenus construct was able to quantify the expression 
levels of Barl, but not the Barl -induced degradation rate of a- 
factor in the extracellular medium, which is why we preferred the 
use of the Fusl-GFP data. 

Third, the a-factor secretion rate could be calculated by 
parameter optimization from images of mixed haploid yeast 
cultures containing MATa Fusl-GFP cells and MATa mCherry 
cells and the information we obtained for the Barl activity. Due to 
the lack of evidence for an extracellular protease activity in MA Ta 
cells, we neglected potential differences in a-factor induction of a- 
factor. In practice, the induction of a-factor secretion should be 
nearly homogeneous on a single image, but may vary for different 
images. We obtained secretion rates from significant fits (F-test p< 
0.05) of 9 images (see Text SI). The mean fitted value of the a- 
factor secretion rate of 865 molecules per second and MATa cell is 
in good agreement with recent experimental measurements (550 
molecules per second and cell measured as basal secretion with 2.5 
- 4-fold maximal induction) [21,22]. The fully parameterized 
model could now be used to efficiently calculate the entire a-factor 
distribution on arbitrary images. We validated the obtained model 
and parameters by predicting the Fusl-GFP fluorescence on an 
image of a larger mixed yeast population not used for model fitting 
(F-test p<2.2e-6). The power of this combined imaging and 
modeling approach is illustrated in Figures 2 and 3 and in Movie 
SI, which shows a mixed haploid population during growth and 
mating next to the simulation of the distribution of a-factor on a 
computational grid generated directly from the corresponding 
microscopic image. 



Global Barl activity limits the range of the a-factor signal 
and optimizes its information content 

We observed large differences in the estimated local a-factor 
concentrations between wild type cell populations and cell 
populations with a barlA background. Dense wild type cell 
populations showed a strongly localized a-factor distribution at 
sites of high MATa cell density, with a-factor concentration 
quickly declining with distance. Consequently, MATa cells far 
away from a cluster of MATa cells experienced significantly lower 
local a-factor concentrations than close-by cells, and hence were 
often non-permissive for induction of the pheromone response 
(Figure 3). Populations with barl A background showed an almost 
uniform distribution of very high pheromone concentrations, 
resulting in global pathway activation as evidenced by high Fus 1 - 
GFP expression. Nevertheless, the global (over-) activation led to 
reduced mating events. 

We wanted to see whether this behavior arises in general and 
independendy of the exact spatial composition of the culture. 
Thus, we performed a computational study using randomly 
generated cell populations mimicking the ones observed micro- 
scopically with varying cell densities (Figure 4). Each virtual 
population was simulated both with wild type Barl secretion and 
in barl A background. We tracked key parameters such as the 
average a-factor concentration, the pheromone gradients per- 
ceived by the individual MATa cells (calculated as the average 
difference in a-factor concentration a cell would sense at its shmoo 
tip and the opposing cell site), and the maximum information 
content of the a-factor distribution. Information (or Shannon 
entropy) quantifies the "surprisal" of a specific a-factor concen- 
tration, i.e. how likely an observed a-factor concentration is given 
an overall a-factor distribution obtained for many cell populations. 

Virtual wild type populations exhibited a strong gain for the 
information content of the a-factor distribution as the popula- 
tion size increases (Figure 4A). This was accompanied by 
increasing a-factor gradients across MATa cells (Figure 4B). In 
contrast, populations not secreting Barl showed information 
contents close to zero (Figure 4A) as well as insignificant 
pheromone gradients (Figure 4B), both independently of 
population density. We noted that the overall pheromone 
concentration remained within a range of up to 20 nM in wild 
type, but in the mutant linearly increased with population 
density (see Figure S6 in Text SI). This observation indicates 
that the gradients and, thus, the reachability of nearby mating 
partners can only be detected faithfully in cell populations 
secreting Barl, particularly in high cell densities. 

Additionally, we simulated various scenarios where a high- 
density subpopulation was placed next to a low-density subpop- 
ulation (Figure 4C-D). Here, the wild type is capable of limiting 
the a-factor distribution to the corresponding subpopulation, 
leaving the low-density subpopulation unaffected by the high local 
a-factor concentration of the high-density subpopulation. The 
same also holds true for random cell distributions (see Figure S7 in 
Text SI). Again, this behavior was not observed in the absence of 
Barl, showing that Barl activity restricts the distribution of a- 
factor. Hence, only subpopulations with high local cell densities 
and small intercellular distances, as required for successful mating, 
were exposed to a-factor concentrations permissive for mating. 

We could validate this model prediction concerning the 
dependency of mating success on culture density experimentally 
by incubating mixed MATa Rpl9A-GFP/M4Ta mCherry popu- 
lations in cell densities varying from 0.5-10 million cells/ml for 
various time-points up to 5 h. We incubated the cells in Petri 
dishes of 36 mm diameter to ensure that the experimental density 
in the resulting cell layer is in agreement with the simulated density 
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Figure 2. Time course of pheromone distribution. Shown are overlay images of Fus1-GFP and mCherry (2A-2C) and the simulated pheromone 
distributions (2D-2F) at three distinct time points. Yellow arrows indicate mating cells after 2.5 hours. Note that the occurrence of mating events 
corresponds well with areas of high pheromone concentrations and gradients. The extracellular a-factor distribution is indicated on the color scales 
bar. Note that a-factor is recorded in a 3-dimensional volume and projected on the 2-dimensional plane related to the image plane. This results in 
units molecules per am 2 or, alternatively, nM-um, as indicated. The full time course is provided in Supporting Movie SI. 
doi:1 0.1 371 /journal.pcbi.1 003690.g002 



(density and distance calculations can be found in Text SI). 
Subsequent cell counting during bead-normalized flow cytometry 
allowed us to quantify the absolute number of each cell type along 
with the number of diploid cells in each sample (Figure 5). Here we 
also observed a strong dependency of diploid formation on cell 
density where efficient mating was predominantly observed in cell 
concentrations higher than 5 million cells/ml, which coincided 
with cell distances permissive for mating (Figure 5A). At the same 
time we could observe a pronounced growth phase taking place in 
parallel with diploid formation (Figure 5B), giving support to the 
finding that there was indeed a variety of different phenotypes 
(mating and growing cells) occurring at the same time in the same 
culture. Figure 5C visualizes the cell densities used in the 
measurements for comparison and Figure S7 in Text SI shows 
the simulated levels of a-factor under these conditions, with 
representation of a-factor distribution and Barl activity distribu- 
tion for one selected cell density. 

Baii enables population growth coordinated with high 
mating efficiency and its absence leads to global cell 
cycle arrest upon pheromone stimulation 

Our results suggested that Barl acts by restricting the activity of 
a-factor to sites, where successful mating is possible, leaving the 
remaining areas free for continued growth. 

In order to test the validity of this prediction and its dependency 
on Barl we observed mating between AIATa cells (here marked 
with Rpl9A-GFP) and AIATa cells (marked with mCherry) and 
quantified their growth rates with FACS analysis (Figure 6) and 



OD measurements in wild type and in barl A populations during 
incubation (Figure 7). 

Mating rates were quantified by the rate of diploid 
formation. To measure the diploid formation rate, we used 
flow cytometry for mixed populations of AIATa. and MATa 
(where we took care to not destroy extracellular gradients 
before the actual measurement) to quantify the fractions of 
MATa., MATa and MATa/ a diploids over time for a fixed 
number of cell counts (Figure 7A). We found no difference in 
the rate of diploid formation between wild type and barl A 
cultures before completion of the first cell cycle (< 1 20 min). 
This observation is in agreement with our results that positive 
effects on the perceived pheromone gradients require higher 
cell densities (Figure 4A,B). However, after passing the first cell 
cycle, the relative fraction of diploids is clearly larger in the 
wild type cultures than in the mutant, consistent with the 
general view that Barl activity helps to reveal the position of 
mating partners [3,13,14] 

Looking at population growth during mating, we found strong 
differences between wild type and barl A cultures (Figure 7B). For 
barl A cultures, the global activation of the pheromone response in 
effectively all MATa cells of the population led to an almost 
complete loss of population growth (Figure 7B,C). This also caused 
a characteristic population phenotype with many pheromone- 
stimulated AIATa cells being significandy larger than normal 
MATa cells and showing multiple mating projections (Figure 7E,F). 
This phenotype was never encountered in unperturbed wild type 
mixtures of MATa and AIATa cells, but could be induced by 
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Figure 3. Microscopic images (left) and derived computational domains and o-factor distributions (right) for BAR1 wild type (top) 

and /foa/'7z/(bottom). The microscopic images are an overlay of the bright field, mCherry (for MATaj and GFP channels (for Fus1 expression in MATa). 
Individual images are given in Figure S1 in Text SI. Computed images show MATam red and MATa in white. The extracellular a-factor distribution is 
as indicated on the color scales. Note different scales for a-factor for wild type (top) and bar! A (bottom). 
doi:10.1371/journal.pcbi.1003690.g003 



swirling them rapidly to inhibit cell fusion (Figure S8 in Text SI). 
Thus, this phenotype appears associated with induction of 
pheromone response in vivo under conditions where a cell cycle 
arrest has been induced but successful mating is inhibited. 

Wild type cultures exhibited significant growth on the popula- 
tion level despite the higher rate of diploid formation and a normal 
phenotype of MATa cells (Figure 7B,D-E). 

The effect of Barl secretion on haploid growth rates was 
even more prominent when looking at the MATa./ MATa ratio 
in the population (Figure 7C). There is no known secretion of 
an extracellular protease described for MATa cells equivalent 
to Barl. Co-cultured wild type MATa cells strongly outperform 
MATa cells in growth during mating to an extent that within 
5 hours MATa is the predominant haploid cell type in the 
population. This cannot be observed in barl J background 
where the MATa/ MATa ratio remains constant, presumably 
because both haploid cell types are equally inhibited in growth. 
In summary, secretion of Barl enables a high mating rate on a 
population level, but also strongly optimizes the population 
growth rate by avoiding unnecessary cell cycle arrest when 
mating is improbable (Figure 7F). 

In order to further validate our findings about the role of 
Barl in the mating process, we mixed MATa cells with 
different ratios of wild type and barlA MATa cells. We 
measured the amount of haploid and diploid cells after 4 h of 
incubation by FACS analysis (Figure 8). For labeling of MATa. 
cells we again used mCherry, whereas MATa BAR1 wild type 



cells were labeled with Rpl9a-GFP and barl A MATa cells with 
Rpl9a-TagBFP2. For equal ratios of wild type cells of both 
mating types at time 0 h we obtained a diploid fraction of 14% 
at time 4 h (leftmost columns). However, MATa cells assumed 
about 23% while MATa cells reached more than 62% of the 
total population, again supporting the observation that a part 
of the MATa population engages in mating and other cells 
continue to grow. 

For equal ratios of wild type MATa and barl A MATa at 
start, we obtained only 3.1% diploids and roughly equal ratios 
of the haploid cells after 4 h. This is in agreement with the 
view that essentially all cells stop growing and start to prepare 
for mating. Microscopic imaging confirms that all cells are 
shmooing under this condition (Figure 8B). When mixing 50% 
of MATa with different ratios of wild type barl A MATa cells, 
we obtain diploid cells of both types at ratios as could roughly 
be expected from the mixes with either wild type or mutant 
MATa. Strikingly, for small ratios (5%) of wild type MATa cells 
at 0 h, we see that barl A MATa mate more frequently (4.5% of 
total population) than in pure mutant mixes (3.1%), indicating 
that they profit from the Barl secreted by wild type MATa 
cells. Figure 8C shows a microscopic image of a mixture of 5% 
wild type MATa and 45% barl A MATa cells together with 
50% MATa at the beginning. Here, the barlA MATa cells 
exhibit clearly lower levels of shmooing compared to the 50%/ 
50% mix in Figure 8B. A few successful mating events leading 
to diploids are indicated by white arrows. 
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Figure 4. Virtual cell populations were randomly generated in order to track the influence of population density on the (/-factor 
distribution. (A) Maximum information content of the a-factor distribution as calculated by entropy (see Text S1) depending on population density. 
(B) Average a-factor gradients (relative front/back difference, see Text S1) for individual cells in the populations. (C and D) Calculated a-factor 
distributions for subpopulations of different densities in wild type and barl A. 
doi:10.1371/journal.pcbi.1003690.g004 



Again, the experiments mixing wild type cells with cells not 
secreting Barl - the barlA or cheater cells - confirm the role of 
Barl which is even supportive for those cells not actively secreting 
it, but experiencing its effect on a-factor levels and gradients. 

Discussion 

The combination of spatial modeling and a series of in vivo 
experiments employing images of mating cells allowed us to 
quantitatively describe the distribution of the pheromone a-factor 
in the intercellular space and the role of the protease Barl in 
shaping the pheromone pattern. W e found that secretion of Bar 1 is 
a highly cooperative mechanism. Haploid AIATa cells secrete 
individually small quantities of Barl molecules. These few 
molecules are quickly distributed across the population by 
diffusion and generate Barl activity that strongly influences the 
diluted distribution of a-factor. For example, a low cell density of 
about 10 6 cells/ml induces a degradation rate of a-factor in the 
range of 10 ' molecules per second which corresponds to almost 
no degradation of a-factor, however, high cell densities can create 
substantial degradation by the Barl activity where single a-factor 
molecules are cleaved within a second after secretion (compare 
Text SI). Thus, individual A'lATa cells need to secrete only very 
small numbers of Barl protease. The global concentration of Barl 
appears to be fine-tuned to ensure a highly informative a-factor 
distribution. This effect is based on the interplay of the secretion of 
a-factor, its diffusion, the activation of Bar 1 transcription and a- 
factor degradation by Barl. Thus, we observed steep gradients in 



the distribution of a-factor only in high densities of AIATa. cells 
where the joint degradation by Barl limits a-factor diffusion 
sufficiently. Our conclusions differ from what has been shown in 
previous publications where the effect of Barl is demonstrated 
either completely theoretical for one or a few cells [13,14,16] or 
experimentally for an artificial setup in a microfluidic device 
without mating partner [15]. First, we observed little to no effect of 
Barl in a setting with only few cells since the Barl activity is 
insufficient to degrade a-factor before it diffuses a large distance. 
Our results also indicate that Barl rather acts on the level of 
subpopulations than on a few individual cells. Second, we also did 
not observe an effect reported earlier that Barl secretion leads to 
"self-avoidance" [15]. Third, we modeled and experimentally 
quantified the induction of Barl due to stimulation with a-factor 
(compared to constant Bar 1 levels as in [1 3] or constant secretion 
as in [14,15]). 

Remarkably, the simple circuit of pheromone secretion and 
degradation by jointly secreted Barl is able to produce highly 
dynamic behavior in yeast populations. In general, the a-factor 
concentration profile recovers the regions where mating has a 
high chance of success, but quickly drops at all other regions, 
thus allowing cells to continue growth when there is only little 
chance for successful mating (see Figure 7E, F). This is further 
regulated by a stimulated Barl production, which depends on 
the extracellular a-factor concentration. Stimulation of Barl 
secretion might be a strategy to adapt the zone of influence of 
a-factor to varying cell densities and numbers of AIATa cells in 
the population. Cells will react to high a-factor concentrations 
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Figure 5. Density dependence of cell growth and diploid formation. (A) Fraction of diploids in mixed populations of different cell densities 
measured over time (each data point comprises between 5.000 and 260.000 counted cells). Indicated is the starting condition at time 0 h. (B) 
Absolute cell numbers of MATa, MATa and diploid cells at 4 h. Note the stronger increase in cell numbers of MATa compared to MATa. (C) Examples 
of intercellular distances for populations with various cell densities. The shown densities correspond to the ones used in the experiments. Density 
calculations are provided in Text SI, simulated distributions of a-factor and Bar1 activity are given in Figure S7 in Text S1. 
doi:1 0.1 371 /journal.pcbi.1 003690.g005 



with strong secretion of Barl, which culminates in a steady 
state permissive for efficient mating (also see Text SI and 
Movie SI). 

A lack of Barl in mixed haploid populations has crucial 
influence on the population phenotype since it leaves many cells 
in a prolonged cell cycle arrest in G 1 phase along with activation 
of the pheromone response. Moreover, the temporally extended 
stimulation leads to larger cells with many mating projections. 
This result is further supported by the observation that MATa. 
cells strongly outperform MATa cells in growth, an effect 
depending on Barl. Thus, Barl is beneficial for growth as well 
as diploid formation because it enables continued growth for 
large parts of the population, but it also provides an enhanced 
ability to interpret the extracellular pheromone signal at sites 
where many cells cluster into locally dense subpopulations. While 
the gradient-enhancing effect of Barl has been reported before, 
we additionally connect it to the requirement of high local cell 
densities [13]. This indicates that the gradient enhancing ability 
of Barl has evolved to take place selectively in dense yeast 
populations and is not a treat of individual yeast cells in large 
volumes. 



Our observations suggest that the yeast populations segregated 
into a spatial pattern with localized regions of diploid formation 
and other regions of continued growth. As a consequence, the 
entire population is divided into two different work programs. 
Both of those programs are performed in parallel and Barl is 
sufficient to induce this separation by forming a locally varying a- 
factor pattern. Depending on the overall cell density, the locations 
reserved for mating show high local concentrations of a-factor, 
whereas the locations reserved for growth show negligible 
concentrations (compare Figures 5 and 7). However, one may 
wonder how this cooperation arose during evolution since cheater 
MATa cells not producing Barl can also easily exploit it. A 
possible explanation lies within the spatial structure of the 
population, since the fitness benefit conferred by Barl is locally 
restricted. Haploids can only arise from a small set of initial spores, 
which form subpopulations by budding and mating type switching. 
As a consequence, cells profiting from the collective secretion of 
Barl are likely to be genetically related. Under these circum- 
stances there is indeed evidence that cooperation can be conserved 
during evolution [23,24]. Our experiments show that cheater cells 
not producing Barl can indeed profit from the Barl secreted by 
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Figure 6. Flow cytometry has been used to quantify the fractions of diploids and MATa or /W/47>/haploids in the yeast cultures. MATa 
and MATa carried constitutively expressed GFP and mCherry constructs, respectively (Rpl9A-GFP and mCherry induced from the TDH3 promoter). We 
observed a switch from many cells carrying either of the two constructs at the initial time point (A, C) to a large subpopulation carrying both 
constructs (B,D) after 5 hours. Density is indicated by colors from blue to red, two biological replicates of 10.000 cells each. 
doi:10.1371/journal.pcbi.1003690.g006 



non-cheater wild type MATa. cells in a mixed population (compare 
Figures 7 and 8). However, the non-cheater cells consistently 
outperform the cheater cells in mating as well as growth, indicating 
that the non-cheater cells maintain an advantage even in a mixed 
population of cheaters and non-cheaters. 

The analyzed regulatory circuit of combined pheromone and 
protease secretion is not only observed in Saccharomyces cerevisiae, but 
is also found in other fungi [25-27]. Furthermore, a similar 
mechanism is known in Dictyostelium discoideum secreting phospho- 
diesterase (PDE) during detection of cAMP [28,29] . Taking this 
into consideration, the described mechanism might be a general 



strategy to separate a cell population into subpopulations with 
different transcriptional programs. 

Our methodology of quantifying the distribution of extracellular 
morphogens in the absence of direct measurement also has 
potential applications in other problems of cellular communication 
and pattern formation. A reduction from the computationally very 
expensive 3D problem (especially for parameter estimation) to an 
integrated 2D problem is feasible for any cells that sediment to the 
bottom of the containing volume under non-agitated conditions. 
However, this computational tool could be used to model the 
behavior of any culture in a non-moving liquid film such as on the 
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Figure 7. The two haploid strains and formed diploids were tracked using flow cytometry of cultures initially containing mixtures 
of MATa Rpl9A-GFP and MATa mCherry cells. (A) The population fraction of diploids in wild type and bar! A tracked over time (compare Figure 
S1 in Text SI, n = 20.000). (B) The population fraction of the two haploid cell types in wild type and barlA cultures tracked over time (n = 20.000). (C) 
Average cell counts for the mixed populations in wild type and barlA (n = 6). (D) Average cell diameters for the mixed populations in wild type and 
bar! A (n = 6). (E) Phenotypes of mixed cultures initially containing MATa Rpl9A-GFP and MATa mCherry after 26 hours. Images are overlays of the 
bright field, mCherry and GFP channels. (F) Illustration of the proposed mechanism for the local coordination of growth and signaling by Bar1. 
doi:10.1371/journal.pcbi.1003690.g007 



surface of fruits or any controlled fermentation such as wine or 
beer production where the liquid is kept still for some time. This 
makes the method applicable for clinical research as well since 
biofilm formation involving quorum sensing is a major complica- 
tion when fighting bacterial infections [30,31]. Here, it might be 
helpful to quantify the distribution of quorum signals in order to 
find possible ways to optimally disrupt the system. 

Materials and Methods 

Used strains and constructs 

Wild type MA Ta reporter strains used in this study are Fusl- 
GFP and Rpl9A-GFP. They are based on BY4741 [MATa. his3Al 
leu2A0 metl5 JO um3A0) and part of the yeast GFP collection [32]. 
The MATa. reporter strain expressing mCherry under control of 
the TDH3 promoter (MATa canlA STE2pr-SpHIS5 lyplA::STE3pr- 
LEU2 kis3Jl leu2A0 ura3A0 metl5A hoA0::TDH3pr-mChmy- 
MATMX4) was a friendly gift of Alexander DeLuna [33]. As 
mutant MATa. reporter strains we used two different strains: barlA 
Fusl-GFP and barlA Rpl9a-TagBFP2. The first mutant was 
created by deletion of the BAR1 gene in the Fusl-GFP strain 
mentioned above, the second was cloned by tagging of the Rpl9a 
Gene with Tag-BFP2 in the BY4741 barlA strain. Rpl9a tagging 
was used because of the high expression level and also since it is 
not known to be involved in the mating process. Yeast strains were 
cultivated at 30°C in synthetic medium (0.17% yeast nitrogen base 
without amino acids, 0,5% ammonium sulfate, 2% glucose, 
55 mg/1 adenine, 55 mg/1 L-thyrosine, 55 mg/1 uracil, 20 mg/1 
L-arginine, 10 mg/1 L-histidine, 60 mg/1 L-isoleucine, 60 mg/1 L- 
leucine, 40 mg/1 L-lysine, 10 mg/1 L-methionine, 60 mg/1 L- 
phenylalanine, 50 mg/1 L-threonine and 40 mg/1 L-tryptophane). 
BAR1 deletions in MATa. reporter strains were inserted by 



homologous integration of a URA3 cassette in the BAR1 locus 
(barlA0::URA3). PCR amplification of the URA3 cassette from 
plasmid template pESC-Ura (Stratagene) was done by sequential 
amplification with the primer pairs 1 / 2 and 3/4 shown in Table 1 . 

This was followed by transformation and selection on agar 
plates with synthetic medium lacking uracil. Verification of the 
BAR1 deletion was done with a physiological assay based on 
growth inhibition by a-factor pheromone [34]. 

The barlA Rpl9a-TagBFP2 reporter strain was cloned by PCR 
amplification of a TagBFP2 loxP-Ura3-loxP transformation 
cassette with primer pairs 5 and 6 from Table 2. As PCR 
template we used vector EKP232. EKP 232 was cloned by ligation 
of TagBFP2 into PstI site of pUG72 [35]. 

The qVenus expressing strain under control of the Barl 
promoter (Bar 1 pr-q Venus) was cloned by using the plasmid pSP 
34 from Serge Pelet [20]. Barl promoter region [—500 bp] was 
amplified from genomic DNA (strain BY4741) by PCR using 
primer pair 9 and 1 0 as well as the first 51 bp of the BAR1 gene 
using primer pair 7 and 8 shown in Table 3. PCR product of 
promoter and gene were mixed and used as template for a fusion 
PCR with a Pmll restriction site between promoter and gene. 
Included in forward and reverse primer were restriction sites for 
SacI and PstI respectively. The fusion PCR product was ligated 
into SacI/PstI site of pSP 34 resulting in plasmid EKP252. 
Plasmid EKP252 was linearized using Pmll restriction enzyme and 
used for transformation and homologous integration into the 
BAR1 locus under preservation of the BAR1 gene. Positive clones 
were selected in minimal medium lacking leucine, a-factor 
induced qVenus expression was controlled microscopically, and 
wild type Bar 1 activity was verified by a physiological assay based 
on growth inhibition by a-factor pheromone [32] and comparison 
with Barl wild type and deletion strains. 
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50 62.2 - - 25 30.6 37.5 46.9 12.5 15 5 5.5 45 56.3 - - 100 99.8 - - 
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Figure 8. Do wild type and mutant cells influence each others' mating success? (A) Mating in mixed cultures initially containing MATa cells 
(marked with mCherry) and varying fractions of wild type MATa cells (marked with Rpl9A-GFP) and MATa barl A cells (marked with RPI9a-TagBFP2) 
after 4 h. (B) Confocal merged image of a mixed culture initially containing 50% MATa and 50% MATa barl A cells. (C) Confocal merged image of a 
mixed culture containing 50% MATa cells, 45% MATa barl A cells and 5% wild type MATa cells. Some mating events are marked with white arrows. 
Each percentage represents the mean of two technical replicates with 10.000 counted cells each. Note that colors in the microscopic images and in 
the table (A) concur. 
doi:1 0.1 371 /journal.pcbi.1 003690.g008 



Confocal microscopy and data analysis 

Microscopic images were acquired with an inverted FluoView 
1000 microscope (Olympus, Tokio, Japan) equipped with a 60 x 
(1.2 N.A) water-immersion objective and a climate chamber 



(Tokai Hit, Japan). GFP was excited with a 488 nm argon laser 
and mCherry with a 559 nm laser diode. Fluorescence emission 
was detected in the range 500-545 nm and 570-670 nm, 
respectively. The Bar lpr-q Venus construct was excited with 



Table 1. 


PCR primer for BAR1 deletion. 






Name 


Sequence 


Use 


Primer 1 


5'-GAAGGGTCATATAATGTCGCGCGTTTCGGTGATG-3' 


PCR 1 


Primer 2 


5'-CTCCAGATTTCTTAGTTTTGCTGGCCGC-3' 


PCR 1 


Primer 3 


5'-GGTTCGTATCGCCTAAAATCATACCAAAATAAAAAGAGT GTCTAGAAGGGTCATATAATG-3' 


PCR 2 


Primer 4 


5'-GACTATATATTTGATATTTATATGCTATAAAGAAATTGTA CTCC AG ATTTCTTA-3 ' 


PCR 2 


doi:1 0.1 371/joumal.pcbi.l 003690.t001 
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Table 2. PCR primer for Rpl9a -TagBFP2 tagging. 



Name Sequence 

Primer 5 5'-GGTATCTACGTTTCTCACAAGGG"nTTATTACTGAAGATTT AGGAGCAGGTGCTGG-3' 

Primer 6 5'-CTGCTACTTTAAAGAAAATGTCACAAAATCAAATAAAAAGC GCATAGGCCACTAGTGGATCTG-3 ' 



doi:1 0.1 371 /joumal.pcbi.1 003690.t002 



515 nm and detected between 530 nm and 630 nm, for Rpl9a- 
TagBFP2 we used 405 nm excitation and as detection range 425- 
475 nm. For mating experiments, MATa Fusl-GFP wild type and 
barlA reporter strains as well as the MATa reporter strain 
(TDH3pr-mCherry), were cultivated to mid logarithmic phase and 
mixed equally. Mating was followed over indicated time periods 
microscopically while microscopic samples were kept in cultivation 
medium at 30°C. 

Image acquisition for a-factor calibration curves was done with 
synchronized cultures of Fusl-GFP wild type and barlA. Cultures 
were synchronized in Gl phase by elutriation with a Beckman 
Coulter JE-5.0 elutriation system. Synchronized cells were 
incubated with a-factor pheromone for 3 hours at 30°C. 
Afterwards cells were spinned down on the surface of a glass 
bottom dish (MatTek Corporation, Ashland, US) by centrifugation 
at lOOx g using self-built accessories. For Fusl-GFP wild type, a- 
factor pheromone concentrations in the range between 0 uM- 
100 U.M were used and for barlA we used 0.1 nM - 1 uM. Mean 
fluorescence intensity of Fus 1 -GFP was analyzed as described in 
the Computational Techniques (see Text SI). For validating the 
employment of Fusl-GFP as proxy for the mating response 
pathway we used a strain expressing qVenus under control of the 
Barl promoter in BAR1 wild type background. Non-synchronized 
cells were incubated with a-factor pheromone as described for 
Fus 1 -GFP BAR! wild type cells and analyzed in the same way. A 
comparison of the results is shown in Figure S5 in Text SI. 

Growth curves and FACS analysis 

Growth of equally mixed MA 7a and MA Ta reporter strains, as 
well as a haploid control strain was analyzed by measuring optical 
density at 600 nm with a Photometer (Eppendorf Bio Photometer 
plus) and in parallel by analysis of cell number and cell size 
distribution with a cell counter (Casy Counter TTC, Scharfe 
System). Yeast cells were incubated in a water bath at 30°C 
without shaking. In time steps of 15 min samples were removed 
from the water bath, vortexed, appropriately diluted, and analyzed 
in duplicate. To quantify the amount of haploid MATa. and MATa 
cells, of diploid cells or of cells within the mating process, we 
measured fluorescence intensities for GFP and mCherry of 10.000 
living cells of each sample by FACS analysis taking advantage of 
the fluorescence of MATa Rpl9A-GFP and MATa mCherry in a 
BD FACS Ariall cell sorter (Becton Dickinson, Franklin Lakes, 



NJ), equipped with a 488 nm and a 561 nm laser with filter sets for 
GFP (525/50 BP, 505LP) and for mCherry (610/20BP, 600LB). 
Cultures were incubated in a water bath at 30°C without shaking. 
In 20 min time steps, duplicate samples were removed from the 
water bath, mixed vigorously, diluted in PBS and FACS analyzed. 
Gates for MATa, MATa and diploids were set by hand identifying 
the cell types as shown in Figure 6. 

Analyzing the influence of the average cell distance on 
the mating process 

As proof of the model prediction we performed a mating 
experiment with different cell densities. A'LATa TDH3pr-mCherry 
and A'lATa RPL9a-GFP BAR1 wild type cells were grown in SD 
medium to mid log phase. Cells were diluted in SD medium and 
cell numbers were adjusted to 10-106 cells/ml by measuring the 
cell number with a CasyTTC cell counter. MATa and MATa cells 
were mixed 1 : 1 and diluted in SD medium in following 
concentrations: 10-106, 5-106, 2.5-106, 1-106, 0.5-106 cells/ml. 
2 ml aliquots of the diluted cultures were incubated at 30°C in 
Petri dishes with a diameter of 36 mm (Falcon), in order to get an 
average monolayer of cells after sedimentation. In time steps of 
1 5 min one Petri dish of each cell dilution was removed from the 
incubator, cells were re-suspended in the medium by intensive 
pipetting and 400 u.1 of the cultures were mixed with 400 |il PBS 
supplemented with CaliBRITE APC Beads (BD Biosciences 
#340487). The samples were analyzed by FACS. APC Beads 
were recorded with 640 nm excitation and 670/41BP filter, 
MATa Rpl9a-GFP and MATa mCherry reporter strains as 
mentioned above. APC beads were gated and used as internal 
standard. In each sample the number of cells corresponding to a 
fixed number of 90000 APC beads was analyzed, giving not only 
the relative amount of haploids and mating events but also the 
growth behavior of the components of the mixed culture (results of 
the experiment are shown in Figure 5). 

Influence of varying amounts of MATa barl A cheater 
cells in mixed cultures with BAR1 wild type cells 

To analyze the influence of barlA cheater cells in mating 
mixtures we used MATa Rpl9a-TagBFP2 barl A reporter strain, 
together with the already introduced MATa and MATa reporter 
strains. The three strains were grown in SD media to mid 
logarithmic growth phase and diluted in SD media to 1-10 7 cells/ 



Table 3. PCR primer for pBarl-qVenus reporter strain cloning. 




Name 


Sequence 


Use 


Primer 7 


5'- GGAGCTCGCGCGAAACTCGCCAA -3' 


BAR1 gene 


Primer 8 


5'- CATCGACACGTGTCTAGAAGGGTC -3' 


BAR1 gene 


Primer 9 


5'- GACACGTGTCGATGAGTCCTTAAG -3' 


promoter 


Primer 10 


5'- CCTGCAGTATATGACCCTTCTAG -3' 


promoter 


doi:10.1371/jou 


rnal.pcbi.1003690.t003 
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ml. Several mating-with-cheaters-mixtures were prepared as 
shown in Table 4. 

Computational methods 

The images were analyzed with CelllD [36] to extract mating 
type, fluorescence activity, as well as position, size and shape of the 
cells. These data were transferred to a computational domain (see 
Figures S2 and S3 in Text SI for details). From this computational 
domain a triangular mesh was generated using Gmsh [3 7] that can 
be used by various RD toolboxes. Here, we used the open source 
Toolbox DUNE [38] to solve the stationary as well as time- 
dependent equations by a finite element method with high 
accuracy [39,40]. In Supporting Text S2 we systematically 
compare 2D and 3D simulations to account for the fact that 
images are taken in 2D, while diffusion and mating happen in 3D 
(see Figures S9, S10, Sll in Text S2). We found that the 
maximum difference between 2D and 3D in the simulated a-factor 
distribution was below 5%. Therefore, for the parameter fit the 2D 
solution was used, which resulted in a major speed-up. 

Supporting Information 

Movie SI Dynamic simulation of an experimental 
mixed haploid population. In the movie generated from 
confocal images to the left, images are overlays of the bright field 
channel, GFP channel (AIATa Fusl-GFP) and mCherry channel 
(AIATa mCherry). The simulated a-factor distribution is shown on 
the right. In the simulation red indicates MA Ta cells and white 
AIATa cells. Each second in the movie corresponds to 5 minutes in 
real time. Due to the delay in the pheromone dependent induction 
of Bar 1, a-factor will initially accumulate to high concentrations, 
which is counteracted by a strong degradation of a-factor caused 
by Barl induction after 30 minutes. Finally, the a-factor level rises 
again due to the loss of Barl induction to its steady state value. It 
should be noted that shmoo formation seems to take place shordy 
after the final stable gradient has been formed, indicating a good 
timing between regulation of the extracellular a-factor distribution 
and mating. The movie requires the Xvid codec (http://www. 
xvid.org). 
(AVI) 

Text SI Computational techniques. Detailed description of 
the image based quantification of pattern formation. Computa- 
tional methods and equations are described in detail as well as the 
employed tools. Additional microscopic images are shown in 
Figures SI and S8. The computational domain and generated 
meshes are explained in Figure S2 and S3. The calibration curves 
referring to Figure 1 can be found in Figures S4 and S5. For 
experimental data and simulations showing the interdependence 
of mating events and cell density see Figures S6, S7 and Table S2. 
Model assumptions and parameters from recent works are listed in 
Table SI. 
(PDF) 

Text S2 Comparison between 2D and 3D model formu- 
lations. Computational quantification of the relation between the 
2D and 3D model for several arrangements of AIATa and AIATa 
cells are shown in Figures S9 to Sll. The transformation of 
important quantities is listed in Table S3 and computational 
results are shown in Table S4. 
(PDF) 
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